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Flower Color Polymorphism in Rhododendron cyanocarpum 
( Ericaceae) , an Endangered Alpine Species 
Endemic to NW Yunnan, China 
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Abstract; Rhododendron cyanocarpum is a narrow endemic species with pink and white floral color. In the present 
study, to investigate the significance of petal color morphs, we examined color morph frequencies, petal color reflec- 
tance and other associated floral characters, effective pollinators, visitation frequencies, and fruit production in the 
field. In all surveyed known populations, plants with pink color morph dominanted and comprised 77% — 100% of in- 
dividuals. Two peaks at 430 nm and 650 nm were found in the petal color reflectance of pink flowers, and only one 
peak at 430 nm was found in the reflectance spectrum of white flowers. In addition, color morphs were also associat- 
ed with colors of style and stigma, lengths of corolla, calyx and pedicel, the closest distance between stigma and sta- 
men, but not with style length and nectar production. Moreover, higher visitation frequency of their shared pollina- 
tors (bumblebees) and fruit production were observed of pink flowers than white flowers. Despite a briefly temporal 
and spacial study, we suggest that color morph frequencies, visit frequencies of bumblebees and fruit production, all 
favor to be stablilizing selection for the pink color morph. 
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Flowers are detected and discriminated by pol- 
linators according to specific signals such as color, 
size, shape, odor and other characteristics. Of these, 
flower color is considered to be associated with polli- 
nator attraction. Therefore the maintenance of floral 
color polymorphism within natural populations is an 
interesting evolutionary puzzle (Frey et al., 2011; 
Rausher, 2008). It is well established that pollina- 
tors are agents of directional selection on floral color 
(e. g. Stanton et al., 1989; Campbell et al., 1997; 
Emms and Arnold, 2000; Wesselingh and Arnold, 
2000; Jones and Reithel, 2001; Streisfeld and 
Kohn , 2005; Frey et al., 2011). Hence, adaptation 
to different pollinators was thought to drive diver- 
gence in flower color and morphology (e. g. Frey et 
al., 2011). Otherwise, if floral color morphs sharing 
the same pollinators, then the color morph most at- 
tractive or frequent should be favored, thereby lead- 
ing to stabilized flower color variation if selection is 
consistent (e. g. Waser and Price, 1981; Smithson 
and Macnair, 1996, 1997; 
Rausher, 2000; Tremblay and Ackerman, 2003 ). 


Moreover, it has been well documented that bumble- 


Subramaniam and 


bees showed a preference for a common floral color 
by frequency-dependent selection when foraging on 
two colors of artificial flowers (Real, 1990; Smith- 
son and Macnair, 1996). 

Compared with interspecific variations, intraspe- 
cific variations of flower color are relatively uncom- 
mon (Matsumura et al., 2006). To explain this dis- 
connect, several authors have proposed alternative , 
but not mutually exclusive, hypotheses where selec- 
tion through pollinators ( direct selection on color ) 
or other environmental factors (indirect selection on 
color) maintains the color polymorphism (e. g. Jones, 
1996; Gigord et al., 2001; Medel et al., 2003; 
Tremblay and Ackerman, 2007; Frey et al., 2011). 

On the basis of field observation, we found pet- 
al color polymorphism in Rhododendron cyanocarpum , 
an alpine species with narrowly distribution and en- 
demic to NW Yunnan, China (Fig. 1A, C). In the 


present study, we investigated color morph frequen- 


cies across all known populations, compared floral 
characters and flower color reflectance , effective pol- 
linators and their visitation frequencies, and fruit 
production among treatments, with the aim of better 
understanding the maintenance of floral color pat- 


terns in R. cyanocarpum. 


1 Materials and methods 
1.1 Plant species and study sites 

R. cyanocarpum is an endangered evergreen shrub 
belonging to Subsect. Thomsonii Sleumer in Subge- 
nus. Hymenanthes (Chamberlain, 1982; Wang and 
Xie, 2004). It is endemic to mountain slopes above 
3 000 m in the Cangshan Mountains around Dali, 
NW Yunnan (Chamberlain, 1982; Wu, 1986). It 
has umbelliform or racemose inflorescences with 5-9 
flowers. Flowers exhibited color polymorphism and 
petal color (pink or white) was consistent within in- 
dividuals (Fig. 1A, C). The flowering period is a- 
bout late March to May. According to three consecu- 
tive years’ observation, flowering time was not diver- 
gent for both pink and white flowering individuals. 
For a single flower, both of pink and white color mor- 
phs flowered approximately 8-10 days. The flowers 
have characteristics of apically porose anthers and ad- 
hesive pollen. The stigma surface usually began to 
become wet at the same time as corolla opening, and 
remained sticky throughout anther dehiscence. 
1.2 Frequencies of petal color morphs among 
populations 

We examined all known five populations of R. cy- 
anocarpum covering both east and west slopes of the 
Cangshan Mountains in April, 2009. In Huadianba 
(HDB) and Ganchaiging (GCQ), we randomly chose 
20 plots with 15 mx 15 m and carefully examined the 
color morphs of all flowering individuals. Due to re- 
stricted distributions, the other three populations, i. e. 
Xiaohuadian (XHD), Guogaishan ( GGS) and Dian- 
shita (DST) , 12 plots with 15 mx 15 m were chosen 
and the color morphs of all flowering individuals 
were carefully examined in each plot (Table S1). 


Then, we calculated proportion of each color morph 
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among total examined individuals per plot. 
1.3 Petal color analysis 

To assess light reflection patterns of both pink 
and white petals at different wavelengths, we ob- 
tained spectral data from petals using a S2000 minia- 
ture fiber optic spectrometer with a PX-2 pulsed xen- 
on lamp (Ocean Optics, Dunedin, FL, USA). All 
measurements were carried out in the range from 250 
nm to 750 nm, using 0.30 nm increments. Differ- 
ences in the reflectance pattern have been found in 
certain plants ( Frey et al., 2011), but preliminary 
testing for changes in the reflectance pattern along 
the same petal did not show any differences, and 
therefore only one measurement per petal was taken. 
In total 20 petals from different flowers, obtained 


from ten individuals for each color morph were ana- 


Reflectance/% 


Reflectance/% 





Fig. 1 


lyzed to assess their potential differences. 
1.4 Floral characters of color morphs 

To determine whether color morph was associat- 
ed with other floral traits, nine floral characters were 
observed at HDB, i.e. style color, stigma color, 
flower number per inflorescence, corolla length, 
style length, calyx length, pedicel length, nectar 
volume and the closest distance between stigma and 
stamen. Flowers were randomly selected from 4 pink 
and 4 white individuals for above mentioned meas- 
urements. To obtain accurate volume of nectar, we 
selected flower buds from 4 pink and 4 white individ- 
uals and bagged to exclude visitors and prevent evap- 
oration (Neg and Corlett, 2000) , and then measured 


it directly from calibrated capillary tube for a volume 


reading when they were newly open (Table 2). 
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1.5 Visit frequencies of the effective pollinators 
to different color morphs 

Effective pollinators have been well documented 
and bumblebees were the mainly effective pollinators 
to R. cyanocarpum (Ma et al., 2010; Ma et al., 2014). 
To investigate visit frequencies of pollinators to each 
color morph, we observed pollinator activity in a 5 m 
x 5m quadrate that only contained one pink flowering 
individual and one white flowering individual, in 
April, 2008. To make the same floral display, each 
plant was left with 300 flowers and we recorded the 
number of pollinators’ visit to the selected individuals 
(from 10:00 to 14:00) in four successively sunny 
days. We then calculated the number of visits to se- 
lected individuals per hour as visit frequency. For sta- 
tistical need, four plots were conducted to newly con- 
structed quadrats under similar weather conditions 
(Sun et al., 2005). In 2009, we used the same 
method and floral display size, and three replications 
were conducted in four successively sunny days. 
1.6 Pollination treatments 

To examine reproductive success under narual 
condition and whether reproductive barrier might oc- 
cur between the two floral color morphs, six treat- 
ments were manipulated from each of 12 inflores- 
cences in 4 pink and 4 white flowering individuals at 
HDB. Totally 76, 78, 67 and 92 flowers, prior to 
anther dehiscence, were selected, then flowers were 
covered with nylon nets after four pollination treat- 
ments: (1) white x white, (2) white ¢x pink @, 
(3) pink d x white and (4) pink x pink, with 
pollen donor from different plants at least 10 m away. 
Other two treatments, 67 white flowers and 75 pink 
flowers in each of 12 inflorescences were randomly 
selected from 4 pink and 4 white individuals as con- 
trol (5), (6). Fruits produced by these flowers 
were counted in October when they were full-sized 
and after fruit abortion had occurred (Ma et al., 
2010). Reproductive success was assessed on the 
basis of fruit set. 
1.7 Data analysis 


All data were examined for normal distribution 


with a one-sample K-S test (Ma et al., 2012). For 
frequencies of each color morph in examined popula- 
tions, data were arcsine-transformed before compari- 
son to normalize the data. One-way anova analysis 
was used to compare frequencies of color morphs, 
floral characters, visit frequencies and fruit sets a- 
mong treatments. Two-way anova was used to assess 
the effects of population and petal color on propor- 
tions of each color morph. All means and standard 
deviations were calculated for all quantitative varia- 
bles using SPSS 11.5 for Windows (SPSS, Chica- 
go, IL, USA). 


2 Results 
2.1 Frequencies of color morphs among populations 
Among the observed 1 711 flowering plants, the 
white flowering individuals only accounted for 7% in 
total (Table S1). Pink and white flowering plants 
were asymmetrically mixed in four populations. The 
white morph plants accounted for 13%, 4%, 3%, 
11%, and pink morph plants for 77% , 96% , 97% , 
89% in HDB, XHD, GCQ and DST, respectively 
(Fig. 2). In GGS, we did not find white color morph 
and the pink color morph hence accounted for 100%. 
Proportion was significant difference between color 
morphs in each of five examined populations (Table 
1). For the same color morph, observed proportions 


also varied among populations (Table 1) , with the 


—@— pink flowered individuals 


—O— white flowered individuals 


Proportion of each color morph 





HDB GGS XHD  GCQ DST 


Fig. 2 Proportions (mean + s. e.) of each color morph 


in the examined five populations 
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lowest proportion of white color morph in GGS and 


highest in HDB (Fig. 2). 


Table 1 Two-way anova table for the effect of population and petal 


color on proportions of each color morph in R. cyanocarpum 





Source of variation d. f. MS F 
Population 4 0. 168 10. 412“ 
Petal color 1 56. 209 3490. 235 ** 
Population x petal color 4 0. 548 34. 01 *™ 
Error 142 0. 016 

Total 152 


“ Significant effect ( P <0. 01) 


2.2 Petal color reflectance 

The reflectance spectrum of petals showed ex- 
tremely low variation between plants ( gray area in 
Fig. 1B, D), all clearly showing a marked peak in 
the reflectance spectrum at 430 nm. However, an- 
other peak reflectance spectrum was also examined 
at 650 nm only in pink flowers. 
2.3 Floral characters of color morphs 


Pink color morph had red style and red stigma, 


but style was lighter green and stigma was pink in 
white color morph. In the observed flower numbers 
per inflorescence, style length and nectar produc- 
tion, no difference existed between two color mor- 
phs. Corolla, calyx and pedicel in white color morph 
were significantly longer than in pink color morph. 
However, for the distance between stigma and sta- 
men, white color morph was significantly shorter 
than pink color morph (see Table 2 for details). 
2.4 Visit frequencies of the effective pollinators 
to different color morphs 

The visit frequencies of bumblebees were 1. 02 
+0. 16 (mean +s. d.) and 0. 73 + 0. 07 in 2008 and 
2009 to pink flowers, and 0.53 + 0. 18 and 0. 50 + 
0. 03 to white flowers in 2008 and 2009. The differ- 
ence is significant in both years (Fig. 3). 
2.5 Pollination treatments 

Fruit set in pink flowers (61%) was significantly 
higher than white flowers (74%) under natural con- 


dition (Fig. 4). For other four treatments, fruit sets 


Table 2 Floral characters of the two color morphs, with the same letter indicating no difference 





. Flower Corolla Style Calyx Pedicel Closest Nectar 
Color Style Stigma ; 2 
i f l number! length length length length distance~ volume 
atterns color color 
P (n=40) (n=20) (n=20) (n=20) (n=20) (n=20) (n=22) 
Pink Red Red 6. 00+0. 99° 4.58+0.55° 3.96+0.42* 0.52+0.09° 1.27+0.10° 0.81+0.08° 0.076+0. 018° 


White Lighter-green Pink 


6.05+0. 88° 4.95+0. 18° 3.96+0.14° 0.68+0.09° 1.81+0. 14° 0.74+0.06° 0.086+0. 017° 


1. Flower number per inflorescence; 2. The closest distance between stigma and stamen 
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Fig. 3 Visit frequencies (mean + SE) of bumblebees 
in each color morph in 2008 and 2009, with the 


same letter indicating no difference 
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Fig.4 Fruit sets (mean + SE) among the six treatments (1, white 
x white; 2, white ĝ x pink 2 ; 3, pink 3 x white 2 ; 4, pink x 
pink; 5 and 6, natural conditions of white and pink flowers) , 


with the same letter indicating no difference 
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(all>87% ) were higher than natural condition treat- 
ments, but are not significantly different within these 


hand cross pollination treatments (Fig. 4). 


3 Discussion 

Few reports have been involved in intraspecific 
variations of flower color mixed in the same popula- 
tion but see Aquilegia caerulea ( Miller, 1981), 
Dactylorhiza sambucina ( Gigord et al., 2001; Pel- 
legrino et al., 2005), Limonium wrightii ( Mat- 
sumura et al., 2006) and Lepanthes rupestris ( Trem- 
blay and Ackerman, 2007). R. cyanocarpum had 
pink and white flowers, and the proportions of two 
color morphs varied among populations, nevertheless 
high proportion of pink color morphs was consistent 
in all surveyed populations. In GGS, no white color 
morph was found, which can be explained mainly for 
such two reasons. First, considering the low propor- 
tion of white color morph in other four observed pop- 
ulations, it seemed that GGS was only composed of 
pink flowering individuals. Second, white color mor- 
ph existed but was not flowering yet during our ob- 
servation, hence we could have missed the white 
color morph. In observed populations, these two col- 
or morph individuals gathered in the same area, and 
once a white flowering individual was found, often 
with many pink flowering individuals closely sur- 
rounding. So, environmental factors such as pH, 
water availability might not contribute to color differ- 
entiation. 

Floral color variation was often associated with 
other floral characters ( Frey, 2007). Our results in- 
deed proved such traits as colors of style and stigma, 
length of corolla, calyx and pedicel. For the closest 
distance between stigma and stamen, autogamy can 
occur easier in white color morph than in red one, 
which will need to be confirmed in the future studies. 
Abundant nectar generally attracted more pollinators 
(see Mitchell, 2004, for a review). No significant 
difference of nectar production showed flower color 
may be the main reason for visit frequencies of bum- 


blebees. 


Floral color offered a visual cue that stimulated 
pollinator sensory systems and that selectively attrac- 
ted certain types of pollinators ( Stebbins, 1974; 
Melendez-Ackerman and Campbell, 1998; Hodges 
et al., 2002). Our results showed both color morphs 
shared the same pollinators and floral color variation 
was not associated with pollinators’ divergency. For 
the same floral display, pink color morph is capable 
of attracting more visit frequencies of pollinators , 
thereby sucessfully leading to higher fruit set under 
natural condition. It should be pointed out that bird 
might be another pollinator in the late flowering of 
R. cyanocarpum on the basis of personal information 
from local people in the study area. If so, the loss of 
petal color reflectance peak which are attractive to 
bird by the peak reflectance between 600 and 700 
nm in white flowers might indicate the loss of bird 
pollination potential ( Shrestha et al., 2013). 

Divergency of floral color may due to adaptive 
selection or non-adaptive mechanisms, such as ran- 
dom mutation ( Epperson and Clegg, 1987; Brad- 
shaw and Schemske, 2003), genetic drift ( Wilson 
and Thomson, 1996) and indirect response to selec- 
tion on pleiotropically related floral traits (Tremblay 
and Ackerman, 2007; Coberly and Rausher, 2008). 
For R. cyanocarpum, we do not understand which 
mechanism contributed to the formation of such floral 
color morphs. However, once floral polymorphism 
existed, which are frequently assumed to be main- 
tained via natrual selection (Tremblay and Ackerman, 
2007). In the present study, effective pollinators are 
both bumblebees, hence hypothesis that pollinator 
divergency maintains floral polymorphism can be re- 
jected. However, even initial visit frequencies of 
both two color morphs are the same, higher frequen- 
cy of pink color morph (92.9%) favored to be fre- 
quency-dependent selection by bumblebees, as a 
reault of learning (Smithson and Macnair, 1996). 
In addition, fruit set also favored to pink color mor- 
ph under natural condition. Therefore , reasons men- 
tioned above, all favor to be stablilizing selection for 


the pink color morph. However, selection can be 
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complex in space and time, other floral characters, 
such as the closest distance between stigma and sta- 
men, may partially facilitate reproduction by autoga- 
my in the future, thereby maitaining white color 


morph within and among populations. 
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Table S1 The investigated populations, mean and stand deviations (SD) of 
flowering individuals of the two color morphs 
Pink White 
Population Location Plot number 

Mean SD n! Mean SD n 
HDB N25°52’, E99°59' 20 19.35 5-17 387 3.3 1.81 66 
GGS N25°51’, E100°02" 12 10. 92 5.12 131 0 0 0 
XHD N25°51’, E100°02’ 12 27. 67 21.72 1589 0.9 1. 06 9 
GCQ N25°52’, E99°58" 20 29. 05 12. 96 581 1.2 1.19 24 
DST N25°40’, E100°06' 12 16. 17 5:31 158 1. 92 1.51 23 
Total 76 20. 63 10. 06 2846 1.46 1.31 122 


1. Number of investigated plants 


